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ABSTRACT 

An algorithm has been developed to solve the following problem: given 
the characteristics of the engine to be used to make afinitebum maneuver 
and given the desired orbit, when must the engine be ignited and what must 
be the orientation of the thrust vector so as to obtain the desired orbit? 


The desired orbit is characterized by classical elements and functions of 
these elements whereas the control parameters are characterized by the 
time to initiate the maneuver and three direction cosines which locate 
the thrust vector. 


The algorithm has been built with a Monte Carlo capability whereby 
samples are taken from the distribution of errors associated with the 
estimate of the state and from the distribution of errors associated with 
the engine to be used to make the maneuver. 
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GLOSSARY OF SYMBOLS 


{x, y, z} — inertial coordinates of the spacecraft 

fj. — gravitational parameter of the central body 
m — mass flow rate of the engine 
V — exhaust velocity of the engine 

e 

m — mass of the spacecraft 

tp — time to start the engine 

t^ — time of engine shutdown (terminal time) 

0 — vector whose components are the parameters being targeted to 
(terminal conditions) 

e - imit vector whose components locate the direction of thrust 
X(t) - the spacecraft state vector at time t 
a — semi-major axis 
e — eccentricity 
'i. — inclination 

RCA — radius of closest approach 
Vp — orbital velocity of spacecraft at RCA 
a — pitch angle 
/3 — yaw angle 

r* PA0E BLANK 
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AN ALGORITHM FOR TARGETING 
FINITE BURN MANEUVERS 


INTRODUCTION 

During orbiting and interplanetary missions it is usually necessary to perform 
maneuvers so as to attain prescribed boundary conditions. Such boundary condi- 
tions can be given in terms of a particular state to be reached or, more frequently, 
in terms of an orbit or trajectory which the spacecraft must be on at the termi- 
nation of the maneuver. For example, during a planetary orbit insertion maneuver 
the spacecraft is initially on a hyperbolic trajectory; it is required that after 
the maneuver the spacecraft be in an elliptical orbit about the planet and that this 
elliptic orbit have prescribed characteristics. 

In most instances the targeting of a spacecraft during a maneuver is accomplished 
using the impulsive algorithm which does not involve any numerical integration. 
However, as advanced missions become more complex, the impulsive algorithm 
may not be sufficient to calculate fuel requirements for a given maneuver. A 
recent mission which underscores this point is that of the Mariner Mars 1971 
spacecraft where the orbit insertion maneuver lasts for roughly 15 minutes. 

According to Robbins [iJ the impulsive algorithm can fail to produce the same 
results as a finite burn algorithm for two reasons: one is the gravity gradient 
effect and the other is the effect due to non-constant thrust vector orientation 
during a maneuver. The gravity gradient effect is the contribution, to the motion 
of the spacecraft during a maneuver, of the time and position dependence of the 
gravitational acceleration. 

Consequently, algorithms which simulate the motion of a spacecraft during a 
maneuver are needed. This need can be satisfied in two ways depending on the 
mission: 

(i) integrate the linearized equations of motion through the burn 

(ii) integrate the nonlinear equations of motion through the burn 

In both procedures the problem of how to treat the thrust vector arises; here, 
there are three options: 
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(iii) keep the thrust vector orientation fixed throughout the burn interval 

(iv) allow the thrust vector to pitch at a fixed rate, or, 

(v) allow the thrust vector to have three degrees of freedom 

This report will present some results obtained from our development of a program 
incorporating (ii) and (iii). 

The problem of reaching prescribed boundary conditions is a matter of determin- 
ing the correct thrust vector orientation and the best time or place to begin the 
maneuver. The problem is statistical in nature because the actual trajectory 
of the spacecraft is never completely known; the filtering of tracking data yields 
an estimate of the actual location of the spacecraft and associated with this 
estimate is a covariance matrix which yields information about the quality of the 
estimate. 

The thrust vector will be determined by using the best estimate of the trajectory 
and a guess at the thrust vector orientation as the initial conditions in the non- 
linear differential equations describing the motion of a thrusting spacecraft. An 
algorithm which minimizes a given performance index is then used to iterate on 
the initial guess of the thrust vector orientation. 

The Monte Carlo option of the algorithm is carried out in the standard way. 
Explicitly the covariance matrix associated with an estimate of the state of the 
spacecraft is assumed to be available and furthermore it is assumed that the 
standard deviations associated with the thrust magnitude error and two pointing 
errors are known. The algorithm first targets the spacecraft to the constraint 
parameters; this yields an initial time to begin the maneuver and the direction 
of the thrust vector. Next the algorithm diagonalizes the state covariance matrix 
and then samples from this matrix and the diagonal covariance matrix associated 
with the engine. The errors obtained from this sampling procedure are then 
added to the spacecraft state and to the model which is being used to simulate 
powered flight. 

It must be emphasized, however, that during the sampling procedure, the thrust 
vector orientation obtained during the targeting phase is associated with each 
member of the sample. This means that no member of the sample is targeted 
to the constraint parameters. In developing the algorithm in this way, the time 
consuming iteration procedure is used only once, when solving for the time to 
begin the maneuver and the thrust vector orientation. 
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The rationale for this procedure is twofold; (i) the engine to be used for making 
the maneuver is assumed to have been sized so that fuel loading optimization is 
not needed and, (ii) with (i) in view, any maneuver is always performed using 
the best estimate of the state of the spacecraft as a reference so that it becomes 
an esoteric study to determine the effect of varying the initial orientation of the 
thrust vector for each member of the sample. 
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CONTROL PARAMETERS 


The equations governing the motion of a thrusting vehicle are well known to be 


X = - ^ X + 1 cos dj cos 6 

r3 m(t) 


m V 


m 

y = — — y + cos \p sin 8 

J.3 m(t) 


fj. ^ . I 

z - - z + s 1 n 0 

r3 m(t) 


m V 


where x, y, z are inertial Cartesian coordinates, is the gravitational parameter 
of the control body, m is the mass flow rate of the engine, is the exhaust 
velocity and m is the mass of the vehicle. The angles and 8 , shown in Figure I 
serve to define the thrust vector orientation. 



Figure I 
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The magnitude of the radius vector of the spacecraft with respect to the origin 
of the coordinate system is given by 


= x"' + y + 2'' 


The mass flow rate is allowed to vary in the program; however, since the engine 
has been sized the burn time is known. Consequently there are four parameters, 
called control parameters, which characterize the maneuver and must be solved 
for. Specifically, the time to begin the maneuver is unknown as is the orienta- 
tion of the thrust vector at the beginning of the burn. 

The orientation of the thrust vector is given by the three quantities, = cos cos 6 
E = cos 0 sin e and E = sin i// . The initial orbit is known and consequently 
the state of the spacecraft on this orbit can be associated with the time to begin 
the maneuver. Once the orientation of the thrust vector is determined at the 
beginning of the maneuver it is held fixed throughout the burn interval. 
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CONSTRAINT PARAMETERS 


The constraint conditions are characterized by the classical orbit parameters 
and functions of these parameters; they represent terminal boundary conditions 
which must be satisfied by the vehicle at the end of the burn interval. 

In particular the set of constraint conditions is composed of the following 
parameters: 

semi-major axis 
inclination 
true anomaly 
long, ascending node 
velocity of perigee 
angular momentum 

During any particular run, a subset of these conditions is usually used. 

A direct solution of the problem is not possible since the constraint conditions 
are non-linear functions of the control parameters, that is, the final state 
attained by the spacecraft is quite critically dependent upon where the maneuver 
is initiated and upon the orientation of the thrust vector. 


eccentricity 
arg. of perigee 
radius of perigee 
arrival energy 
time from perigee 
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METHOD OF SOLUTION 


The overall problem falls into the general category of two point boundary value 
problems. The technique used in this program to find a solution is to iterate 
on the initial conditions until the given terminal conditions are satisfied. Jn 
particular, let tj be the terminal time, x(t) be the state at time t and 0 [x(tj)J 
the vector of terminal conditions which must be satisfied. In order to satisfy 
this terminal vector, the initial time t^ and the initial thrust vector orientation, 
e(t ) = (E E , E ), must be corrected and recorrected until the terminal 
conditions ’arrsatisfied. To construct this iterative^rocedure it is realized 
that e(t„) completely determines a solution because x(tQ) is fixed once t^ is 
chosen. Thus x(tf ) depends and this implies that the terminal conditions 

are representable in the form© [ • 

In order to formulate the orbit transfer strategy a performance index is employed 
as a criteria to determine a solution. In this program the index is the weighted 
squared distance, with constraints, in the parameter space between the final 
orbit and the desired orbit; the objective is to find the time t^ and orientation 
e(tQ) which minimizes this performance index. 

The algorithm which performs the minimization is known as MINMAX [2] to 
which the interested reader is referred. A very brief overall description of 
the algorithm is presented below. 

Denote the desired terminal conditions by ©j, [e(to), tf,] and define the vector 
y(t) = [ e(t), t] . Let t . be the ith value of the ignition time and e(t^ j ) the ith 
value of the orientation of the thrust vector. A first order Taylor series expan- 
sion about t^., e(t^,) yields 




0] Cy(toi+i> - y('^oi)^ 


where 


'b© b© 


B© B© B© B© 


t 

BEj 3E2 BE3 Bt 


is an m X 4 matrix with m being the dimension of the vector of terminal condi- 
tions. Thus the original nonlinear problem has been replaced with a linear 
one whose solution is direct; however this solution does not satisfy the nonlinear 
problem. It is for this reason that a sequential procedure must be utilized. 
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The sequential process will now be constructed. A scalar W [y(t^.)l is defined 
as follows: 




where W is a weighting matrix introduced so as to make the components of 
the vector 0 - 0 - (V_ . &)^'y compatible among themselves and superscript 

T denotes the transpose. 

It may happen that the minimization of '1' [y(t^j)l will yield a_Ay which is quite 
large; this can occur when 0^ - ^ is large or when '^y( ® ill-conditioned. 

Large values of Ay can easily take us out of the re^on of linearity resulting 
in wild fluctuations of the residual vector 0^ (t^.) -0(t„i). To prevent this from 
happening a constraint is imposed in the form 

(Ay(t^.))-^W^ Ay(t„.) < S, 

where is a square matrix of scale factors and Sq is a constant chosen in 
such a way that at no point during the iteration is the linearization (the Taylor 
series at each step) violated. 

The performance index may now be defined to take the form 

J[y(t„i)] ='I'[7(t^.)] + ^(Ay(t„j))T Ay(t^.) 

The problem now is to find a set of J^t^. ) i = 1, 2, . . . such that 

and such that the sequence { converges to a point ylt^) which minimizes 

J[y]. 

Minimization of the performance index is accomplished by taking the derivative 
of J with respect to (A^ and equating this derivative to zero. This procedure 
yields 
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= - [v_ 

'■ T(toi) 


0]^ W, (0„ - 0) 


t \W^ Ay(t„,) = 0. 
and solving for Ay(t^.) yields 

Thus 

y(toi+i) " y(toi) + ^y(toi) 
and the next iteration can begin. 


11 



NUMERICAL INTEGRATION 


The nonlinear equations of motion are integrated using FNOL2, [3] . This 
routine uses a fourth order Runge-Kutta or a fourth order Adams-Moulton 
method to solve a system of coupled ordinary differential equations. As currently 
modified, the routine uses a double precision arithmetic throughout. Further- 
more, the truncation error can be held within input bounds by giving the user an 
option of automatically varying the step size. 
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ILLUSTRATIVE EXAMPLE 


A l una r orbit insertion maneuver has been used to illustrate the targeting 
procedure. 

The algorithm was used to target the spacecraft using a set of burn durations 
ranging from 5 seconds to 20 minutes. In each case the initial spacecraft mass 
was taken to be 317.37 kg and the exhaust velocity was taken to be 2.84 km/sec. 
The mass flow rate in each case had to be adjusted before making a computer 
run so that the thrust magnitude would be sufficient to attain a lunar orbit. 

To obtain results the following procedure was implemented: a nominal earth- 
moon trajectory was generated and the spacecraft was targeted off of this tra- 
jectory to the constraint parameters, -b = 88f7, RCA = 2428.88 km. In addition, 
for each burn interval considered, a minimum eccentricity for the lunar orbit 
was sought. For each case this targeting procedure yields a thrust vector 
orientation and time to begin the burn. 

The thrust vector orientation, although computed in terms of inertial coordinates 
as mentioned earlier, is printed out in terms of the geometrically more mean- 
ingful pitch and yaw angles shown in Figure II. 
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The pitch angle is the angle measured from the spacecraft velocity vector V to 
the projection of the thrust vector onto the orbit plane. The yaw angle is the 
angle measured from the projection of the thrust vector onto the orbit plane to 
the thrust vector. 

Results from several cases are shown in Table I. 


Table I 


Burn 

Duration 

(secs) 

Values of 
Constraint 
Parameter 

Values Obtained 
from the 
Algorithm 

Pitch 

Angle 

(Deg) 

Yaw 

Angle 

(Deg) 

No. of 
Itera- 
tions 

5 

-1. = 88?76 
RCA = 2428.88 km 
e = 0 

^ = 88.°76 
RCA = 2428.88 km 
e = .0006 

-178.5 

-.002 

8 

20 

Same as above 

n. = 88.°76 
RCA = 2428.88 km 
e = .0004 

-179.1 

-.002 

7 

60 

Same as above 

A. = 88.°76 
RCA = 2428.83 km 
e = .0009 

+179.3 

-.003 

4 


Although the time to initiate the burn was also one of the control parameters, it 
was found in all cases, the time tg never varied by more than 0.4 seconds from 
the initial guess. The reason for this behavior is that the mass flow rate had 
been adjusted properly beforehand to insure orbit insertion. 

The circular orbit velocity at the point of insertion into lunar orbit is 1.4207 km/sec. 
Although the velocity at perigee was not one of the constraint parameters, this 
velocity was calculated for each case; for the three cases shown in Table I, the 
velocity attained is 1.421 km/sec. 

Confidence in the program has been enhanced by a very favorable comparison 
with a program vhich integrates the linearized equations of motion. The compari- 
sons were carried out for burn durations of 5 seconds and 30 seconds; such 
burn intervals are short enough that the linearization of the equations of motion 
does seriously affect the final nominal (targeted) parameters. Such comparisons 
have shown that the targeted parameters (e, RCA) agree to within less than 
half percent. 
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SUMMARY 


A program has been developed which integrates the equations of motion of a 
thrusting spacecraft in order to optimally attain prescribed constraints. The 
algorithm is optimal in the sense that it minimizes a performance criteria which 
is the weighted squared distance between the final constraints and the desired 
constraints. 

The program has been constructed so that a full error analysis of a maneuver 
can be carried out; this capability is available in the Monte Carlo mode of 
operation. 

There are no restrictions concerning the class of orbital maneuvers the program 
will handle so that orbit to orbit transfers as well as orbit insertion maneuvers 
can be investigated. 
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